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Abstract 

We present in this paper a nonlinear sigma-model analysis of a spin-1 
antiferromagnetic Heisenberg chain in an external commensurate stag- 
gered magnetic field. After rediscussing briefly and extending previous 
results for the staggered magnetization curve, the core of the paper is a 
novel calculation, at the tree level, of the Green functions of the model. 
We obtain precise results for the elementary excitation spectrum and in 
particular for the spin gaps in the transverse and longitudinal channels. 
It is shown that, while the spectral weight in the transverse channel is 
exhausted by a single magnon pole, in the longitudinal one, besides a 
magnon pole a two-magnon continuum appears as well whose weight is 
a stedily increasing function of the applied field, while the weight of the 
magnon decreases correspondingly. The balance between the two is gov- 
erned by a sum rule that is derived and discussed. A detailed comparison 
with the present experimental and numerical (DMRG) status of the art 
as well as with previous analytical approaches is also made. 

PACS numbers: 75.10.Jm, 75.0.Cr, 75.40.Cx 

1 Introduction 

In recent times there has arisen a great experimental Q-Q and theoretical 
[D~@l interest in a new class of magnetic materials of the general composi- 
tion i?2BaNiC>5, with R one of the magnetic rare-earth ions (typically: R =Nd 
or R =Pr). These materials are obtained by substitution from the reference 
compound Y2BaNiOs, a highly one-dimensional compound with negligible in- 
teractions among the spin-1 Ni 2+ linear chains. For this reason Y2BaNiOs 
is generally considered as an almost ideal example of an S — 1 Haldane-gap 
system, with a spin gap of: Ao = 0.41048(2) in units of the AFM intrachain 
exchange coupling || M . 

As the magnetic R + ions order antiferromagnetically below a certain Neel 
temperature Tjv, the i?2BaNiOs's have been modelled, to a first approximation, 
as a set of S = 1 chains with a negligible interchain coupling (as compared with 
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the intrachain one), acted upon by an effective commensurate staggered field 
[|| roughly proportional to the sublattice magnetization of the R 3+ lattice, and 
hence increasing when the temperature is decrease below TV. 

The above experimental scenario has motivated a renewal of theoretical ac- 
tivity on the model, that is nonetheless already quite interesting "per se", of 
an integer-spin AFM Heisenberg chain coupled to an external commensurate 
staggered held, that can be described by the model Hamiltonian: 

H = Y,i JS * ' S *+i + (-!) 4H s ' Si} (1) 

i 

where: Sf = S(S + 1) with S an integer (we set h = 1 henceforth, and take 
5 = 1 for the Ni 2+ chains), J > and H s is the external staggered field in 
appropriate units (see, e.g., Ref. Q for details). 

An extensive DMRG study of the model of Eq.(l) has been performed in 
Ref.[||, where very accurate results were reported for the staggered magnetiza- 
tion curve, the spin gaps, the static correlation functions and the correlation 
lengths in both the longitudinal and transverse (with respect to the direction of 
the field) channels. The authors of Ref. || made instead an analytic study of the 
model beginning with the familiar mapping [llj of the Hamiltonian (1) onto 
a nonlinear sigma-model (NLcrM) . What they discussed very accurately was ac- 
tually a related and somewhat more phcnomenological model in which the strict 
NLcrM constraint is softened, then replacing the original NLcrM with a theory 
of the Ginzburg-Landau type parametrized by an adequate set of adjustable 
parameters (see Ref.[|| for details). 

One of the purposes of the present paper is to investigate carefully what 
are the resulting similarities and/or differences when the NLcrM constraint is 
not softened but enforced consistently at each level of approximation. We will 
report here only results at the tree-level of a loop expansion jl2), i.e. essentially 
at the the mean-field (MFT) level, of the partition function of the model, sup- 
plementing them however with a stability analysis, and deferring a systematic 
evaluation of higher loop corrections, that are considerably more involved, to a 
forthcoming paper p3| . 

Our second purpose is to assess the validity in the present context of the 
so-called single-mode approximation (SMA) that did prove beyond doubt its 
validity previously but in rather different contexts |l4| . We will provide here 
the explicit proof of the fact that indeed the SMA is not applicable to discuss 
the elementary excitation spectrum in the longitudinal channel, a claim that we 
had already put forward some time ago [[7| , without providing however there an 
explicit proof. 

The paper is organized as follows. In Sect. 2 we state the essentials of the 
general formalism and derive the saddle-point approximation in the presence of 
a general external source field. This is needed in order to set up the consistent 
scheme of calculation of the propagators at the mean-field level that is reported 
in Sect. 3. In Sect. 4 we study the analytic structure of the propagators, and 
notably of the longitudinal one, at the physical saddle-point, i.e. when the source 
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field becomes the staggered static field of Eq. (1). The final Sect. 5 is devoted to 
a discussion of our results and to a detailed comparison with previous theoretical 
approaches, as well as to a discussion of some as yet unsolved problems that are 
posed by the experimental scenario that has been outlined at the beginning. A 
useful sum rule for the propagators of the NLcrM field is derived in the Appendix. 

2 Saddle-point approximation for a general 
source field 

Under the Haldane mapping Si « (— lfSni + Ij, nf = 1, where 

represents the slowly-varying local staggered magnetization and L is the lo- 
cal generator of angular momentum, the Zeeman term of Eq.(l) becomes : 
EiC-l) 1 ^ • Si « ■ H s + E l (- 1 ) 1 ! t • H s . In the continuum limit 

the second term becomes a total derivative that can be neglected if we adopt 
periodic boundary conditions on the chain. 

Going then to the continuum limit, integrating out the fluctuation field 1 and 
implementing the NLcrM constraint n 2 (x) = 1 (x = (x, r) with r the Euclidean 
time) with the aid of a Lagrange multiplier A = A(x) , we obtain the partition 
function: 3 =Tr{exp[— /3H]} of the model in the continuum limit as the path- 
integral: 



[On] 



eMSeff) (2) 



2tt 

where the effective action S e ff is given by: 

S ,„./^£ eW - S H,»W-,A W (»'W-1,} (3) 

where: J dx = J dx dr and the Euclidean Lagrangian is given by: 

Ce{x) = 2^ (c2|&n|2 + |9rn|2) (4) 

and the NLcrM mapping predicts: g — 2/S for the coupling constant and: 
c = 2JSa (with a the lattice constant) for the spin-wave velocity. 

Now we promote 3 to a generating functional 3 [J] by replacing S e f / with: 

S[J] = J cix{£ £; (x)-S'J(x)-n(x)-iA(x)(n 2 (2:)-l)} (5) 

and we will set: J = H s only at the end of the calculations. 
Altogether (after an integration by parts): 

S[J] = ^J dxrfx'n(x) • G c 7 1 (x,x')n(x') - S J dxJ(x) • n(x) +i J dxA(x) 

(6) 
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and Go solves, with the appropriate (Matsubara-Bose) boundary conditions the 
equation: 

--(c 2 d 2 +d 2 + 2z.gcA(x))G (x,x') =<5 (2) (x-x') (7) 
Performing now the linear shift: n(x) = n'(a;) + a(x), with: 

a(x) =S J rfx'G (x,x')J(x') (8) 



we obtain: 



S[J] = \j dxdx'n'(x)-Go 1 (x,x')n'(x')+i J rfxA(x) 
-\s 2 j dxdx'J(x)-G (x,x')J(x') 



(9) 



Now we can integrate out the field n', obtaining: 
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exp(-S[A;J]) (10) 



3[J]oc 
where: 

S[A;J] - |Tr{In(Go 1 )}+* J dxA(x) - ^S 2 J dxdx'J(x)-G (x, x')J(x') 

(11) 

We analyze now what are the general features of a saddle-point approxima- 
tion made in the presence of an arbitrary (space-time dependent) source field 
J(x). We will also evaluate here the propagators at the mean-field level. The 
saddle point will be determined by the equation: 



( SS[X;3] 
V *A(x) 



= (12) 

j 



where (..)j means that we (functionally) differentiate while keeping J constant. 
As: 

SG ° x u , f ) = 2*G (x',x)Go(x,x") (13) 
oA(x) 

and: 

*^y-3G (x,x) + 5 2 | dydy'G (y, x)G (x, y')[J(y) • J(y')] - 1 (14) 
we find the saddle-point equation in the form: 

3G (x,x) + S 2 y rfyrfy'G (y,x)G (x,y')[J(y)-J(y')] = l (15) 
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Eq.(15) will determine then a space-time dependent saddle-point that will be a 
functional of J as well: A = A* [x; J]. 

In Mean-Field Theory (MFT) one approximates 3 [J] as: 

3[J] «cxp(-S[A*;J]) (16) 

and hence the connected two-point propagators: 

Gf (x,x') = S> {<T(»»(*)»V)> - (»"W> <»V)» = J^jf^ 

(17) 

will be given by: 

G Q/3 fx x") S3 S2S i X *: J ] (18) 

where now, when taking functional derivatives, one has to consider not only the 
explicit dependence of S on J but also the implicit one through A*. We find 
then: 

8S f f 8S \ SXjj) f SS \ 
But the first term on the r.h.s. vanishes identically at the saddle point, and so: 

jj^ = («^i) A - 4 s2 / <,y[Co(x ' y)+0 ° (y ' x)1J " (y) (20) 

which provides also the mean-field equation for the (staggered) "magnetization" 
induced by the source field J. 

Proceeding one step further we find eventually: 

Gf(x,x') = is 2 [G (x,x') + (x^x')]<5 Q/3 

+S 2 J dyrfy'[Go(x,y)G (y,y') + (x^y')]J a (y) (*^fy)( 21 ) 

This is the general structure of the mean-field propagators for a general 
source field J. 

The functional derivative of A on the r.h.s. of Eq.(21) is determined by: 



Sjaysxjj \5X5X'J J 5J a \SJ a V ^A 



J/ A 



where the second term represents the variation of (SS/5X)j w.r.t. its explicit 
dependence on J. As all the quantities in brackets have to be evaluated at the 
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saddle point, this equation becomes an inhomogeneous linear integral equation 
for SX/SJ a whose kernel is: 

= 6r(x,x')+45 2 G (x',x) J dydy'Go(y,x')G (x,y')[J(y)- J(y')] 

(23) 

where T is the "polarization bubble" : 

r(x,x') = Go(x,x')G (x',x) (24) 
The integral equation for <5A/<5J Q (x) reads then: 



J dyH( X ,y) (*^^y) = -2^ 2 Go(x,x') J rfyG (x, y) J a (y) 



(25) 



In the next Section we will specialize the results obtained here to the physical 
case J = const. = H s , which will determine the physically relevant saddle point. 

3 Results at the physical saddle point 

When J = const. = H s , the associated saddle point will correspond also to 
A = const. Setting then: —2igc\ = const. — c 2 £~ 2 , translational invariance 
will be restored and the following results can be easily derived 0: 

1. The saddle-point condition is: 

3Go(0) + 5 2 H 2 [G^)(0)] 2 = 1 (26) 

where Go(0) = Go(q = 0) (q = (q,fl n = 2im/j3) and Go(q), the Fourier 
transform of Go(x), is given by: 

Go(q) = ni + /{q> + ^) (27) 

whence: Go(0) = g£ 2 /c. Explicitly, at T = 0: 

| in { ac + vTTTAeF} = 1 - (?) 2 H ' e4 (28) 

and the cutoff can disposed of by fitting it to the zero-field gap Ao = c£ -1 
that is know from the DMRG studies. Eq.(28) will determine then the 
field dependence of the correlation length £, and it is clear that £ will 
depend quadratically on the field. 



G 




2. The magnetization is given by: 

m s = S (n s ) = S 2 A(0)H s = ^-^H s (29) 

c 

Comparison with the DMRG data of Ref . shows a slight overestimate of 
the values of the magnetization for small fields, but the agreement becomes 
better and better as the field increases (see Fig.l). 

Note that, in view of this equation, the saddle-point condition can be 
written also as: 



3Go(P) + {-f) =1 (30) 

3. While the transverse susceptibility is given by: \ T = m s/H s = gS 2 ^ 2 /c, 
the longitudinal one is given by the full derivative: \ L = dm s /dH Sl i.e. 
by: 

xW { 1+W0 J_}^{ 1+ M} (31) 

The derivative on the r.h.s. can be obtained by differentiating the saddle- 
point equation w.r.t. H Sl and the explicit expression for \ L is: 



x L = x T 



7 



(clearly exhibiting: \ L < X T always, the two coinciding only when H s = 
0). In the limit A£ > 1, Eq.(32) reduces then to: 

%L = xT l + ^(2gSeH s /c)2 (33) 

4. In the translationally-invariant case, Eq.(25) for 5X/5J a can be rewritten, 
using the equation for the magnetization, as: 

J dyH(x y) (*^7y) - -2m?G (x - x') (34) 



Making then the "Ansatz" : 

. 5A(x) 



(5J Q (x') 



m^X(x-x') (35) 



(Hence: SX/SJ a = in the directions orthogonal to the field) Eq.(34) 
reduces to the following equation for X: 

J dyH (x - y)X(y - x') = -2G (x - x') (36) 

that can be solved by Fourier transforming it, thus yielding: 

X(q) = -2^M (37) 
ff(q) 

where, now: 



ff(q) = 6r(q)+4(-^J G (q) (38) 

5. The longitudinal connected propagator is given by: 

G c L (x-x') = 5 2 G (x-x')-2m2 J d y G (x - y)X{y - x') (39) 
while: 

6. The transverse propagator (that has no nonconnected parts) is simply 
given by: 

Gj(x-x') = S 2 G (x-x') (40) 
Going to Fourier space: 

G^q) = S 2 G (q) (41) 
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which is (see Eq.(27)) a free boson propagator that, when analytically 
continued to the real axis, has simple poles at: uj = is(q) with: 



e(«) = V c V + Ag, A = cC 1 (42) 

Notice also that: 

X T = GT(0) (43) 

Explicitly, in Fourier space, one finds: G^(q) = 5 2 Go(q) — 2TO 2 G (q)A'(q) 
and, with some algebra: 

^(q) = ^G 0(q) -^-^)_ = g* (q? 3f (q) 



3r(q) + 2(m s /S) 2 G (q) 3r(q) + 2(m s /S) 2 G (q) 

(44) 

T(q) is the convolution of two Go's that can be evaluated explicitly as: 

~ 1 f Hk 1 

T(q) - -{gcf / ^coth(/3e(fc)/2) 



2 va ' J 2tt v v " 'e(k)e(k + q) 

e(k + q)+e(k) e(k + q) - e(k) 



ni + {e{k + q)+e{k)Y m + {e{k + q)-e{k)f 



(45) 



The analytic continuation in frequency (f2„ ^ iz) has a branch-cut along the 
entire real axis for all (3 < +oo. The discontinuity across the branch-cut vanishes 
however exponentially with temperature in the range: — 2Ao(q) < Re(z) < 
+2A (<?), where: A (q) = c^/(q/2) 2 + = e{q/2) and, right at T = 0, the 
branch cut extends only from — oo to — 2Ao(g) and from +2Ao((?) to +oo. The 
reason for this is that the second term in curly brackets vanishes exponentially 
for [3 i— > oo. The same term vanishes also, irrespective of temperature, for 
q = 0. For this reason we expect also T(q) to be an essentially positive quantity 
at all temperatures. According also to Eq.(27), this implies that H(q) is a 
positive-definite quantity, and hence that _ff(x — x') is a positive-definite kernel. 

Now, the kernel H determines actually also the stability of the saddle point. 
Indeed, by expanding the action (11) quadratically around the (physical) saddle 
point in terms of: 5A(x) = A(x) — A*, we obtain 

S[X] = S[X*) + i J dxdx'5\(x)H{n - x')^A(x') (46) 

and the positivity of H will then guarantee that the physical saddle point is 
indeed differentially stable, i.e. a local minimum. 

We analyze now the asymptotic behaviour of the longitudinal propagator for 
both large and small values of |q|. 

While one sees immediately that: A(q) ss |q|~ 2 for large |q|, in the same 
limit p|: f (q) w ln(|q| 2 )/|q| 2 , instead. Therefore: <Sf (q) w 5f (q) for large 
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|q|, and the (very) short-distance behaviour of the two propagators is the same. 
This implies: 

limGf'(x) = limG^(x) = S 2 G Q (0) (47) 

Recalling that, in the translationally-invariant case, the full longitudinal prop- 
agator G L is related to the connected one by: G L = + mf, we see that, as 
a consequence: 3G (0) + (m s /S) 2 = S- 2 {2G T (0) + G L (0)}, and that therefore 
the saddle-point condition can be again read simply as one implementing the 
constraint on the average, i.e. as: 

(n 2 > = 1 (48) 
At the opposite end, when q i— ► 0: 



~ = 3^r ( o)G (o L 

3r(0) + 2(m s /^)2G (0) 

that is markedly different from GjT(O), coinciding with the latter only for H s 

0. The same will be true for the small-momentum behaviour of G^(q). The 
long-distance behaviours of the two propagators will be then definitely different 
for H s ^ 0, and so we expect quite different asymptotic behaviours at infinity, 

1. e. quite different correlation lengths ||. 

On top of that, the equation (cfr, Eq.(43)): \ L = Gc (0) provides us also with 
an independent expression for the longitudinal susceptibility, explicitly showing 
how it results from both one- and two-magnon contributions (the latter coming 
from the "polarization bubble" T(q)). 

In explicit terms, we have, at T = 0: 

oo 

f (o) = M / d k4^ = l dk{e + r 2 r« = ^ (50) 



87r J e(fc) 3 4irc J Aire 



i) 



Inserting this expression (together with the known value of Go(0)) into the equa- 
tion for Gg(0) yields back precisely Eq.(33) that had been obtained by letting 
A n oo in the previous equation for the longitudinal susceptibility wherever 
this did not lead to divergent results, which is precisely the attitude that has 
been taken in the present calculation. All this shows that (at least) there are 
no inconsistencies in the MFT approach that has been adopted here. 

In the next Section we will discuss in detail the structure of the analytic 
continuation of the propagators to the complex frequency plane and to the real 
axis. 



10 



4 The analytic structure of the propagators 



At the present level of approximation the transverse propagator is just a free- 
boson propagator. Analytic continuation is straightforward, leading to: 



+ . 



(e(q) — cy q 2 + £~ 2 ). Going to the real axis from above: z = u> + ir], 77 1 — ^ 

ImGf (q, w) = {5(u - e(q)) - 6(u> + e(q))} (52) 

The spectral weight function is then fully exhausted by single poles at lu = ±e(q) 
which is the structure required ^, [lj] for the applicability of the SMA. The 
relation, which is a direct consequence § of the SMA: X t = Sgc/{A T ) 2 , with: 
At = Ao = c£ _1 , is obeyed exactly, at this level of approximation, in the 
transverse channel. 

Let's turn now to the longitudinal propagator, and let's begin by looking 
at the polarization bubble T(q) = T(q, f2 n ). We shall consider for simplicity 
only the T = limit in which the second term on the r.h.s. of Eq.(45) can be 
neglected. Then: 

f (q) = \{ 9 cf [ ^ m I , , n7 ^^^^l^ (53) 



2^ ' J 2tt e(k)e(k + q) tl 2 n + (s(k + q) + e(k)) 2 
that we write for short as: 



where: 



and: 



Note that 



E(k,q) =e(k + q)+e{k) (55) 



E(k,q)=E(-k-q,q) (57) 

and the same will hold true for A(k,q). 

Performing now the analytic continuation and going to the real axis from 
above: 

T(q,w) = T 1 (q,uj)+iT 2 (q,uj) (58) 
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where: 



('"P" standing for Cauchy principal part) and: 

T 2 {q 1 u J ) = - A {gcf J dk ^ - - [5{w - E(k, q)) - 5{u + E(k, q))] (60) 

that is odd in u> (which implies that Ti(q,uj) will be an even function of w) and 
positive for positive u>. 

To evaluate T 2 explicitly it will be enough to consider the positive frequency 
part. Notice that, for any fixed q: mm{E(k,q)} = 2A (q) (A (q) = e(q/2), 
A o (0) = A ), and therefore: T 2 (q,w) = for \ui\ < 2A (q), as we know already. 
Otherwise, it is easy to see graphically that the equation: u) = E(k,q) has two 
solutions at k = k*(q,ui) = fc — q/2 and at k = —k* — q = fc — q/2, where: 
fc - M2cW(cj* - 4A^)Vp - c 2 q 2 ). 

Then we obtain easily, as: e{k)e(k + q)(dE(k,q) / dk) = c 2 [ke(k + q)+ 
(k + q)e(k)}: 

r 2 (g,w) = Sgf^ - 4Ag(g))*snM (6i) 

where: C(q) = \k*e(k* + q) + (k* + q)e(k*)\ and: k* = k*{q,\u\). An ex- 
plicit analytic expression for T 2 can then be written down in general, but it 
is not especially illuminating, although it can be very useful for numerical cal- 
culations. It simplifies greatly for g^O, where we get simply: C(q = 0) = 
(M/2c)VS*=4A* 

The (integrable) square-root singularity at the edges of the branch cuts 
is present at finite q as well, and indeed, for uj 2 > 2Aq((/), wc obtain, to 
leading order: C(q) w a(q)y/u 2 - 4A§(<7) + 0(uj 2 - 4A§(qr)) with: a(q) = 
(A (g)/2cAo){|a;| - g 2 c 2 /2A (g)}. 

By exploiting the parity of I^, Ti will be given then, via dispersion relations, 

by: 

OO 

ri(g,w) = 2 I ^JV 2 {q,uj')V{^-^) (62) 

2A ( 9 ) 

which shows that T\(q, oj) will be strictly positive for \uj\ < 2A (q), a result 
that will prove to be useful shortly. 

Let's consider now the full longitudinal propagator: 

G^,n n) = 9 cs 2 ^ m^m (63) 

(3T(q,n n )/2gc)ini + e*(q)] + (m s /S) 2 
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Performing the analytic continuation, omitting specification of the label q 
and defining: 

G(z) = G^(q,fl n )/gcS 2 , T(z)=3T(q,n n )/2gc, S=(m s /S) 2 , e = s(q) 

(64) 

we are led to study the analytic structure of a function of the form: 

^> " m J% ) + s < 65 > 

(0 < S < 1, e < 2A ). T(z) will be given by: 

r ( *)=/^M m 



with T 2 (u>) having all the properties that have been listed above (odd in u, 
positive for positive w and vanishing for \u\ < 2A (cfr. Eq.(61)), thus producing 
a branch-cut in T(z) for real z = u> and 2A < \u>\ < +oo). We will write: 
T(uj + i0 + ) = Ti(oj) + T 2 (u)) on the real axis. The analytic properties of G(z) 
will be determined in turn by its spectral weight function. 
Going to the real axis we find, for z = u + iff, f] > 0: 

G(w + irj) = Gi(w)+»G 2 (w) (67) 

where, defining: 

A{uj) = (e 2 - lu 2 + ry 2 )ri(w) + 2r/ujT 2 (Lu) + S (68) 



B(uj) = - (e 2 - uj 2 + r, 2 )T 2 {uj) (69) 

Giand G 2 are given by: 

Gl {uj) = ^h + b 2 h (70) 

and: 

G2M = — ^ H+B 2 H — ( 71 ) 

Even for 77 1— » 0, G2 needs not vanish when T 2 does. In particular, let's 
inspect its structure for \ui\ < 2Ao. Sending r\ to inside the T's (and only there) 
we have, for |w| < 2A : A(u) = (e 2 - uj 2 + ?7 2 )ri(w) + 5, B(u) = 2^ x{uj) , 
leading to: 

r = ^-^ + ^ + 5/^) 

U ; (e 2 -w 2 +r ? 2 +5/r 1 (a;)) 2 + Vw 2 V ; 
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and: 



G 2 H 



2rjuj 

(e 2 - oj 2 + r, 2 + S/T^w)) 2 + 4t] 2 uj 2 



From the structure of G 2 (uj) it is clear that: 

limG^u;) = Trsgn(uj)S(f(uj)); f(cj) = oj 2 — e 2 — 



7JI->0 



Remembering that Fi is an even function of ui, we can write: r x 
and f(uj) = will have solutions at oj — ±€l with: 



6 2 + 



ri(e|) 



that reduces (as it should) to e 2 

do; 2 



£ 2 (</) when S 



2oj 



S 



* 0. Moreover: 

dT^LU 2 ) 



According to known formulas, then: 



where: 



lim G2 {oj) 
771 — >o 



7 : 



j—{5(u; 



T 2 (lu 2 ) dto 2 



e L ) - S{lu + e L )} 



_5 dTtju; 2 ) 

T 2 (lu 2 ) diu 2 



(73) 



(74) 



Ti(^ 2 ), 



(75) 



(76) 



(77) 



(78) 



This proves of course that the longitudinal propagator has (in the range we are 
examining) simple poles on the real axis at: u = ±€l, and the prefactor 7 will 
give the reduction of the quasiparticle weight w.r.t. the pure bosonic case. 

As, for small fields, 5 oc H 2 , 7 will approach 1 quadratically in the field 
when H s 1— > 0. A numerical plot of the relative quasiparticle weight 7 at q = 
is presented in Fig. 2. It is a steadily decreasing function of the field, and the 
quadratic regime near H s = is confined to a very narrow region of fields. For 
higher fields there is an intermediate region in which 7 is almost linear in the 
field, and we find numerically that for S = 1 it saturates in the high-field limit 
to: 



lim 7 = lim 7 ~ 0.279 



(79) 



More than 70% of the quasiparticle weight is then lost when the system 
evolves towards saturation. No significant changes are expected for q 7^ 0. 

To complete the analysis we have to investigate the range \u)\ > 2A of (?2(w) 
where T 2 (w) does not vanish. With some long but straightforward algebra we 
find: 



G 2 (io) = 



(e 2 - u 2 ) 2 (T 2 (uj)) 2 + (S+ (e 2 - u; 2 )r 1 ( W )) 2 



(80) 
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Figure 2: The relative quasiparticle weight of the pole of the longitudinal propagator as a 
function of the staggered field and at q = 0. 



for \lo\ > 2Aq. Therefore, in this range of frequencies Gi will vanish as H s t— > 
and we will recover the simple pole structure with the longitudinal and 
transverse propagators becoming equal. The longitudinal pole will survive also 
up to saturation, but with a strongly field-dependent strength. 

That as the field increases the spectral weight that is lost from the pole gets 
transferred to the two-magnon continuum (80) (and viceversa when the field 
decreases) is dictated, e.g., by the sum rule: 

J ^-LuImG 2 (Lo) = 1 (81) 

— oo 

The sum rule (81) is just one of the general sum rules connected with the 
moment expansions of the spectral weight functions that are related to equal- 
time expectation values of multiple commutators and that have been known for 
a long time in many-body theory [[IB], [lTj . A proof of the sum rule adapted to 
the specific context of the NLcrM will be given in the Appendix^. 

The pole at represents the longitudinal magnon. It will be a well defined 
excitation as long as e\(q) < 4A§(<7), which we will prove to be the case. It 
will be higher in energy (as we have proved previously that Ti(q,uj) > for 

< 2Ao(q)) than the two (degenerate) transverse magnons that both have 

1 Note however that while (cfr. Eq.(52)) all moments exist for the transverse propagator, 
when G2 is given by Eqs.(77) and (80) only the first moment will exist and all the others will 
turn out to be divergent. This is just an artifact of the mean-field approximation (see, e.g., 
the discussion of a similar problem in Ref . Wn . 
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energy e, and will become degenerate with the latter when 5 i— » (the limit in 
which the full SO (3) invariance is restored). 

Just as in the transverse case (e = e(q) = \Jc 2 q 2 + A T , At = c^ 1 = A ), 
we can define a longitudinal gap Al via: Al = £l(<? = 0), i.e.: 



Al = A T + 



ri(o,A L ) 



In general: 



1^(0,00 = ^1™ 



du' 



LO' 2 - UJ 2 



2A„(q) 



7T V(^' 2 - 4A2) (uj' 2 +uj 2 + r, 2 ) 2 - 4uj 2 uj' 2 



(82) 



(83) 



As we are assuming here: A^ < 2 At, the limiting procedure becomes trivial 
and we obtain: 



ri(o,A L ) 



du> 



1 



2At 



t V(^ 2 - 4A* ) ^ 2 - A| 



Eq.(82) can be rewritten in dimensionless form as: 

2 



where: 



Al 
A t 



= 1 + 



46 



3gF{A L /A T ) 



F(y) 



dx 



7T \/x 2 - 4 x 2 - y 2 



ltfl<2 



(84) 



(85) 



(86) 



Explicitly: 



— — tan 



27TJ/V1 - y 2 /4 1 2 

and, defining: y = At/ At we obtain, with: g = 2/5: 



-Vl-2/ 2 /4 



i + 



8m2j/5 ^1 - y 2 /4 



Sjl-ftan- 1 [fyrr 774]} 



(87) 



(88) 



The numerical results for the longitudinal gap are reported for S = 1 in 
Fig. 3, where we report also the results for the transverse gap. Here too the 
results are in excellent qualitative agreement with the DMRG results, with some 
small quantitative discrepancies in the low-field regime. We also find that near 
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Figure 3: Results for the spin gaps in the longitudinal and transverse channels: 
-Dots and full line: DMRG results of Ref. [^j and our mean-field results for the transverse gap. 
-Triangles and dot-dashed line: DMRG results of Ref. [Ej and our mean-field results for the 
longitudinal gap. 

-Dashed line: SMA predictions for the longitudinal gap (Ref.ja] and our Ref. |7l). 



H s = the longitudinal gap increases roughly as three times the transverse one 
(actually: lim^ sM0 (^L — A )/(Ar — A ) ~ 3.58). For H s i— > oo instead: 



v z = i 



8ySy/l - y 2 /4 



and, for 5 = 1: 



{l-ftan- 1 [f yi-F74]} 



Ar, 

lim — - ~ 1.855 



(89) 



(90) 



i.e. for high helds Al tends to increase slightly less than twice At- This is in 
agreement with the DMRG results of Ref.||. Notice that, however, the limiting 
form of Eq.(89) tells us immediately that the ratio of the saturation values of 
the gaps will tend exactly to two in the large- 5 limit, and indeed it is not 
difficult to see that: 



lim — — 

i?«H-> 00 At 



a 
52 



(91) 



for S 3> 1, with a a numerical constant of order 0.5. The relative quasiparticle 
weight 7 can also be shown to be of order S~ 2 in the same limit, i.e. it will 
vanish when the magnon poles reaches the edge of the continuum. 
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5 Discussion and conclusions 

We summarize here our results, comparing them at the same time with those 
obtained within other approaches. We will list and discuss only a few relevant 
points, namely: 

i) In zero field the excitation spectrum consists of the well-known degenerate 
triplet of massive Haldane bosons with energy e(q) = c\J q 2 + £ -2 with a gap 
Ao = c£ -1 . For finite fields, instead: 

ii) The staggered magnetization curve (Eq.(29) and Fig.l) turns out to agree 
well with the DMRG results of Ref . || . As can be deduced directly from Eq. (29) 
and as was discussed in more detail in Ref. Q , the low- and high-field behaviours 
of the staggered magnetization are respectively: 



with A a numerical constant (see for more details) for large H s . The value 
obtained in R] of \ T — 23.74/ J is somewhat higher than the DMRG result || 
of x T — 18. 5/ J. We will resume this point shortly below. 

Eq.(93) shows that the staggered magnetization saturates only asymptoti- 
cally in the large-field limit. This is what had to be expected, of course, and is 
a direct consequence of our implementing in a consistent way the NLcrM con- 
straint. The authors of Ref.[[| found instead a magnetization curve that (see 
Fig.l) agrees better with the DMRG data than ours in the initial part (i.e. for 
low fieds), but that disagrees more and more as the field is increased. Even 
worse, their magnetization saturates at a finite value of the staggered fields 
and keeps growing indefinitely beyond saturation afterwards. Full saturation 
implies of course that the fully polarized Neel state should become the exact 
ground state of the Hamiltonian (1), which known to be true only asymptoti- 
cally for H s i— ► oo. Moreover, an cwersaturated magnetization is clearly not an 
admissible result. 

It appears therefore that softening the NLcrM constraint as done in Ref.Q 
can be an admissible procedure in the zero-field limit, that can be also given 
(see Jl8| and references therein) some legitimation in the same limit, but that 
can become more and more dangerous as the field increases, leading eventually 
to rather unphysical results. 

Hi) The transverse channel is saturated by a single well-defined magnon pole 
with energy er(q) = s(q) and a gap A T = A .The strength of the pole has a 
rather weak field dependence and the pole persists up to the highest fields. The 
relation: \t = Sgc/ (At) 2 , which is typical of the SMA, is exactly obeyed in 
this channel at the mean-field level. Both our results and those of Ref. || agree 
quite well with the DMRG results. 



m s « X T H S + 0{H 3 S ) 



(92) 



for H s ss and: 



m 




A 



(93) 
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iv) In the longitudinal channel instead we find also a well defined magnon 
pole with energy cl(q) and a gap > At, but the magnon propagator acquires 
also a (two-magnon) branch-cut for nonzero staggered fields. The strength of 
the pole of the propagator at cl is strongly field-dependent and as the field 
increases it is steadily transferred to the continuum, in agreement with the sum 
rule (81). Near saturation, and for 5 = 1, the strength of the pole is reduced to 
less than 30% of its zero-field value, while it vanishes completely at saturation 
in the large- 5 limit. In the same limit the longitudinal pole disappears into the 
continuum at saturation, while it remains below the continuum (hence a well- 
defined excitation, although with a strongly reduced intensity) for finite 5. The 
relation analogous to case Hi), namely: \l — 5c/c/(Ai) 2 is badly violated in 
this case (cfr. Fig. 3) except in the zero-field limit, and this proves that, due to 
a non negligible contribution from the two-magnon continuum, the SMA is not 
applicable at all in the longitudinal channel in the presence of a nonvanishing 
(staggered) field, as was done in ^ (see the close critical comparison with the 
DMRG results that was made in Ref. ||, and the SMA curve that is reported 
for comparison in Fig. 3). This makes a great difference with the zero- field case, 
where there is quite convincing evidence |L9| that two- and/or multi-magnon 
excitations carry a negligible spectral weight, which turns out to be actually 
exactly zero at the mean-field level. 

In conclusion, we would like to point to two somewhat unsatisfactory aspects 
of the theoretical scenario that has been outlined here and elsewhere [|)-|7|. 

First of all, and as we have already remarked, there is excellent qualita- 
tive and even semiquantitative agreement between our NLerM analysis and the 
DMRG results over the whole range of fields. There remain however some quan- 
titative discrepancies. As to these, we can only stress the fact that the NLoM 
mapping is basically a semiclassical expansion starting from the large-5 limit 
that is then continued to lower values of the spin. Being so, and if it has to 
have any validity at all, it should definitely be able to capture all the essential 
features of the low-energy physics of models such as that of Eq.(l), which we 
believe we have proved to be the case for 5=1. It should then lead definitely 
to qualitative and semiquantitative agreement, but it cannot be expected to 
yield also, and it would be actually a totally unexpected bonus if it did, precise 
quantitative agreement with, e.g., the DMRG results. We believe that, apart 
from other refinements []l3f , not much more than that can be expected from the 
NLcrM mapping. 

In view of what has been just said, however, the quantitative agreement 
should improve for increasing values of the spin. A test of the NLcrM, in con- 
junction with a DMRG analysis, on higher-spin chains can be of some interest 

Coming now to the experimental scenario outlined in the Introduction, while 
essentially all the theoretical models presented so far (including ours) yield quite 
reasonable agreement between theory and experiment as far as the staggered 
magnetization curve and the transverse gap are concerned, recent measurements 
[Q of the longitudinal gap in Nd2BaNiOs point to some discrepancies between 
theory and experiment. Namely, the longitudinal gap is found to survive, with 
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Ai > At, also above the Neel temperature, i.e. also when H s = 0, which 
points to (single-ion) anisotropies that, although could be accomodated very 
easily within the NLcrM approach, thus leading to an additional spin-gap in the 
longitudinal channel, are not explicitly accounted for by the models considered 
so far. Moreover, below Tjy, the longitudinal gap appears to grow more slowly 
than the transverse one, and this remains a so far unresolved puzzle. 

More measurements on compunds of the class Ji^BaNiOs are called for to 
ascertain whether this is a general fact or it is specific to Nd2BaNi05. If the for- 
mer were the case, this would imply that a model involving just a single chain in 
an effective staggered field, interesting in itself and worth of theoretical interest 
in its own, although it can account for most of the experimental observations, 
is not enough to provide a complete description of all the observed low-energy- 
physics features of such compounds, and that more sophisticated models are 
needed. 



6 Appendix. A sum rule for the NLcrM 
propagators 

To cope with the notation adopted in Sect. 4 after Eq.(64), let's consider (setting: 
g — c = S — 1), in real time now, a classical Lagrangian of the form: 

£(x,t) = ~{\d t nf - \d x n\ 2 - V(n)} (94) 

with V(n) an arbitrary interaction term, e.g. one implementing in some appro- 
priate limit [|l5[ |l6| the NLcrM constraint n 2 (x,t) = 1. The canonical momen- 
tum density is defined by: 

dC 
d{d t n) 

and the Hamiltonian density will be given, as usual, by: Tt(x,t) = ir-dtn — C, 
with canonical (equal-time) Poisson brackets (PB's): 

{n a {x,t)^{y,t)} = 5^5{x-y) (96) 

with all the other PB's vanishing. 

Canonical quantization (here too we set: h — 1 and use for simplicity the 
same notation for the classical variables and the corresponding quantum opera- 
tors) is accomplished by replacing the nonvanishing PB's with the commutators: 

[n a (x,t),^(y,t)]=i8^5(x-y) (97) 



20 



Then, we obtain immediately, taking expectation values in the ground state: 

-id t {[n a (x,t)y(0,0)])\ t=0 = -*(K(x,0), 7/(0,0)]) - -S(x) (98) 

As the r.h.s. involves only equal-time commutators, this relation will be true 
irrespective of the form of the interaction term as long as the latter depends 
only on n and not on its time derivative. 

Introducing the Fourier transform A a ^(k,uj) defined via: 

(K(x, t), n"(0, 0)]) = f £ / ^ exp{t[kx - ut)}A^[k, co) (99) 



Eq.(98) will imply the sum rule: 



/ 



< ^wA a P{k,u) = 5 a P (100) 



Let's test it in the case: V(n) = m 2 n 2 , with m fixed, e.g., by the condition: 
(n 2 (ir, t)) = 1. Then the field can be quantized Jl8| asf]: 

n a (x,t) = J 4 ^ fc |afc exp[z(fca; - e(k)t)} + a\ exp[— i{kx - e(k)t)]j (101) 



(e(k) = \/k? + m 2 ), and: 

[n a (x, t), ir^ly, t)} = iS a ^S(x — y) & [a%, of ] = 4ir6 aP e(k)6(k - q) (102) 
One finds then easily: 

/Ah 1 
^ e Mikx)— y[exp(-i E (A0t) - exp(ie(fc)t)] 

(103) 

whence: 

A Q/3 (fc, w) = tf^-ZL^w - £ ( fc )) ~ 5 ( w + £ ( fc ))l ( 104 ) 

E(k) 

and the sum rule can be checked at once. 

To compare with Eq.(81), let's recall first of all that our definition (Eq.(17)) 
of the Euclidean propagators differs by a sign from that that is usually adopted 

. Accordingly, we will define the retarded Green functions of the components 
of the n field as: 

Gf(x,t) = i6(t)([n a (x,t),n f> (Q,Q)]) (105) 



2 Note that, in order to reproduce correctly the commutation relations (97) between the 
fields and the conjugate momenta, the signs_of the arguments of the exponentials had to be 
taken here as the opposites of those of Ref. . 
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and their Fourier transform are then given by: 

G£{k,u) = - / —— \s^ + 106 

J 2n uj — uj' + 10 

Hence: 

lmGf{k^) = ^A^{k,J) (107) 

In view of the standard relationship ^2) between retarded (in real time) 
and causal (in Euclidean time) Green functions, the sum rule that we have 
just derived extends to the Euclidean propagators introduced in the text and 
becomes then the sum rule (81). 
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